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Project Objective 

The objectives of this research project were: 

1) To study the fluid dynamics of sheared particle-liquid suspensions and the impact 
of differential particle-fluid inertia 

2) To develop new techniques for observing suspension particle contact and 
deposition upon solid surfaces. 

Project Description 

Dr. Yoda was supported by the NASA Office of Biological and Physical Research on a 
four-year grant from March 2000 through November 2004 for a ground-based study on 
the fluid dynamics of sheared particle-liquid suspensions and the impact of differential 
particle-fluid inertia on such flows. Such inertial effects can only be observed in 
reduced-gravity environments since they are overwhelmed by buoyancy effects on Earth. 
Moreover, these inertial effects will have a significant impact upon suspension flows in 
microgravity. Suspension dynamics are of importance in a wide variety of advanced life 
systems applications, including water reclamation and dust mitigation in confined 
habitats. 

Project Significance 

This study establishes the experimental methodology and techniques required to carry out 
micro- or reduced-gravity experiments on the impact of differential particle-fluid inertia 
on particle-fluid interactions and suspension flows. These effects can only be observed in 
reduced-gravity environments since they are overwhelmed by buoyancy effects on Earth. 
Moreover, these inertial effects will have a significant impact upon suspension flows in 
microgravity. Suspension dynamics are of importance in a wide variety of advanced life 
systems applications, including water reclamation and dust mitigation in confined 
habitats. 

Project Accomplishments 

A new aqueous low-viscosity (viscosity |i = 5 cP) refractive index- and density-matched 
suspension has been developed consisting of non-colloidal polymethyl methacrylate 
(PMMA) particles suspended in a ternary mixture of water, glycerin and ammonium 
thiocyanate. 1 The maximum transmission through this suspension T (defined as the ratio 
of the illumination power transmitted through 1 cm of suspension to that transmitted 
through 1 cm of the suspending liquid) for a volume fraction <j> = 10% of particles with 
radii a - 50 pm is about 60% at 22.5°C for illumination at 488 nm. Comparisons with an 
organic solvent-based suspension system consisting of PMMA particles suspended in a 
ternary mixture of UCON oil, Triton-X and dibromohexane used for laser-Doppler 
velocimetry studies 2 show that this suspension system has comparable or slightly better 
transmission characteristics. Given its relatively nontoxic nature and low viscosity, this 
suspension system will be useful for a variety of suspension mechanics experiments. 

The dynamics of neutrally buoyant prolate and oblate spheroids immersed in a simply 
sheared liquid were studied using a plane Couette channel as a simple three-dimensional 
model of suspended particle dynamics in a dilute suspension (i.e., in the absence of 
particle-particle interactions). Recent two-dimensional numerical studies using the 



lattice-Rnltzmann technique have suggested that these bodies, when constrained to rotate 
about one of their principal axes, undergo a rotation transition fro 1 1 j a time-periodic to a 
stationaiy state above a critical shear-based Reynolds number and that this transition is 
independent of particle shape except for a sphere in simple shear flow. 3 More recent 
three-dimensional lattice-Boltzmann simulations suggest, however, that the dynamics of 
neutrally buoyant spheroidal particles in (plane) Couette flow — when free to rotate about 
any axis — are more complex, with prolate spheroids undergoing two rotation transitions 
and nutation between these transitions, while oblate spheroids undergo only one rotation 
transition. These new results demonstrate that the classical problem of a sphere in 
simple shear flow is actually a special case that cannot be generalized to other particle 
shapes, and that even slight deviations from sphericity due to surface roughness (for 
example) may have a significant impact on suspension particle dynamics. 

A plane Couette flow facility based upon a previous design 5 , was designed and built to 
study cylinders and spheroids in simply sheared single-phase liquids and suspension 
flows. The facility, which is compatible with the new aqueous, density- and refractive 
index-matched suspension system, consists of an endless belt stretched over two large 
rollers immersed in the working fluid; the belt is then driven by a 90V % hp DC motor. 
The flow coordinate system is defined with x along the flow direction, y along the 
velocity gradient direction, and z along the vorticity direction. The test section half-width 
(i.e., y-dimension) h = 1.27 cm; the test section length) and depth (i.e., x- and z- 
dimensions, respectively) were 20 cm and 16 cm, respectively. The linearity and 
steadiness of the simple shear flow the test section were characterized using both single- 
component laser-Doppler velocimetry (LDV) and particle-image velocimetry (PIV). For 
conditions corresponding to a particle- and shear-based Reynolds number Rea = a l G / v < 
20 (a = particle or spheroid semimajor dimension; G = shear rate; v = fluid kinematic 
viscosity), the linearity of the flow, quantified by the deviation between the average 
measured and the expected x-velocity components normalized by the belt speed £/ b , was 
less than 6%. 



Figure 1 Side view of plane Couette facility for generating simple shear flows for both liquids and 
suspensions. Note the two large rollers flanking a white circular cylinder in the center of the test section. 


Rotation periods T were first measured for neutrally-buoyant circular and elliptical 
cylinders with aspect ratio AR = b/a = 0.75 (a and b - cylinder scmimajor and semiminor 
dimensions, respectively) immersed in a simply sheared suspension using our aqueous 
suspension system at various shear-based Reynolds numbers. Figure 2 plots the average 
rotation period normalized by the shear rate, as a function of Re q (for the suspensions, v 
was calculated using the Krieger model 6 ) at a confinement ratio k = a/h = 0.25. These 
measurements are compared with previous results for single-phase flow. 7 The symbols 
represent different particle volume fractions <j>: <j> = 0 (O), 0.01 (O), 0.03 (A) and 0.06 
(□), with open and filled symbols denoting circular and elliptical cylinder results, 
respectively. The data for cylinders in simply sheared liquids and suspensions are in 
reasonable agreement. 





Figure 2 Log-log plot of normalized rotation rate Q / G vs. shear-based Reynolds number Re G for circular 
and an AR - 0.75 elliptical cylinders in simple shear flow of suspensions at volume fractions up to 6% and 
single-phase liquids. 

For the three-dimensional studies, spheroids were manufactured from three different 
materials for these studies: 

1) DMS Somos 7545 resin (specific gravity SG = 1.20) was fabricated into spheroids 
using stereolithography rapid prototyping; 

2) Polydimethylsiloxane (PDMS) (SG = 1.07) was molded into spheroids using 
molds manufactured with stereolithography rapid prototyping; 

3) Acrylic (polymethylmethacrylate, or PMMA) (SG = 1.17) was milled into 
spheroids using a micro-CNC (computerized numerical control) mill. 

In all cases, the spheroids had an aspect ratio AR = bla = 0.33 and 0.5, and a ~ 0.45 cm, 
giving a confinement ratio in the test section k = a/h = 0.35. A variety of spheroids with 
different SG were developed so measurements could be carried out for both single-phase 
and suspension flows. 






Rotation periods were measured for spheroids constrained to rotate about both their major 
and minor axes. Figure 3 presents rotation period data for die two prolate spheroids with 
AR = 0.33 [top] and 0.5 [bottom] rotating about their minor [left] and major [right] axes. 
The symbols denote different spheroid materials, namely DSM Somos (triangles) and 
PDMS (circles); the open and filled triangles denote different types of belt drive to attain 
different ranges of Re g- No rotation transition was observed for both types of prolate 
spheroids rotating about their major axes, presumably because the transition (if it occurs) 
occurs at a critical Reynolds number Re a beyond the range attainable with the plane 

Couette channel (i.e., Re a > 35). A transition to a stationary state was observed, 

however, for both prolate spheroids rotating about their minor axes, with Re ct =19.2 and 
33.5 for AR = 0.33 and 0.5, respectively. 


Prolate Spheroid (AR = 0.33): Minor Axis Prolate Spheroid (AR = 0.33): Major Axis 
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Figure 3 Normalized rotation rate as a function of shear-based Reynolds number for various prolate 
spheroids. 

Figure 4 presents similar data for the two oblate spheroids with AR = 0.33 [top] and 0.5 
[bottom] rotating about their minor [left] and major [right] axes. No rotation transition 
was observed for the oblate spheroids rotating about their minor axes, again most likely 







because the transition (if it occurs) occurs at a critical Reynolds number Re CT beyond the 
range attainable in these experiments. A transition iu a stationary state was observed, 
however, for both oblate spheroids rotating about their major axes, with Re cr = 10.3 and 
30.5 for AR = 0.33 and 0.5, respectively. 


Oblate Spheroid (AR = 0.33): Minor Axis Oblate Spheroid (AR = 0.33): Major Axis 



Figure 4 Normalized rotation rate as a function of shear-based Reynolds number for various oblate 
spheroids. 

The prolate and oblate spheroid cases that do show a transition closely follow the scaling 
law proposed by Ding and Aidun, 3 with GT [Re a — Re G ] as Rea — > Re cr . Most 
interesting, although these bodies were constrained to rotate about one of their principal 
axes, we observed this transition at Re CI = 33.5 for an AR — 0.5 prolate spheroid rotating 
about its minor axis, in direct contradiction of the lattice-Boltzmann results reported by 
Qi and Luo 4 , who reported that this body would remain rotating about its minor axis for 

R e G < 50. We are currently checking our results, since this discrepancy may be due to 

friction in the spheroid mounting in the experiments. We then plan to write up the 
spheroid rotation period results for publication in Experiments in Fluids. 
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